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^ ' The early steps of light response occur in the outer segment of rod and cone photoreceptor. They 

O ■ involve the hydrolysis of cGMP, a soluble cyclic nucleotide, that gates ionic channels located in the 

■ outer segment membrane. We shall study here the rate by which cGMP is hydrolyzed by activated 

phosphodiesterase (PDE). This process has been characterized experimentally by two different rate 
constants Pd and (3sub'- Pd accounts for the effect of all spontaneously active PDE in the outer 
segment, and Psub characterizes cGMP hydrolysis induced by a single light-activated PDE. So far, 
no attempt has been made to derive the experimental values of /3d and /3sui, from a theoretical model, 
which is the goal of this work. Using a model of diffusion in the confined rod geometry, we derive 
analytical expressions for f3d and Paub by calculating the flux of cGMP molecules to an activated 
PDE site. We obtain the dependency of these rate constants as a function of the outer segment 
Q I geometry, the PDE activation and deactivation rates and the aqueous cGMP diffusion constant. Our 

formulas show good agreement with experimental measurements. Finally, we use our derivation to 
model the time course of the cGMP concentration in a transversally well stirred outer segment. 
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I. INTRODUCTION 



I The modern theory of chemical reactions originates back to Arrhenius [1], who showed in 1889 that the 

OO . backward rate constant of two reactants depends exponentially on the temperature and the activation 

energy barrier. However, the molecular description of the backward rate started with the seminal paper 
. of Kramers in 1940 [2| (see also [1, 3)- The constant fcf, is used to describe the chemical reaction of 

abundant species in solution and the concentration of the product resulting of the interaction of two 
molecules is calculated by using the mass action law. But, The computation also involves the forward 
kf rate and the concentration of the two species. At a molecular level, kf reflects the mean time for one 
of molecule to meet the other by diffusion, and the probability to react upon encounter. For diffusion 
limited chemical reactions, based on the mean time for a uniform concentration of particles inside an 
infinite 3-dimensional space to hit a sphere of radius a, von Smoluchowski obtained in 1914 the first 
, estimate kf = AnaD Q. The Smoluchowski formula was later on extended to the case of a partially 

■ absorbing sphere [1, 0] ■ 

Diffusion plays in many cases a prominent role in the determination of the forward binding rate 
[Tol . [ill . [T^ , and numerous fundamental processes in cellular biology rely on the rate at which diffusing 
molecules hit a small target site: examples are trapping in patchy surfaces [Tsj . receptor dwell time inside 
a synapse and many more. When the number of molecules is not large, the mass action law is not 
sufficient to account for the random nature of the chemical reactions and other approaches are required 
p^ . In addition, in a confined geometry, the Smoluchowski formula does not describe the refine structure 
of the bounded space. For that purpose, the small hole approximation was developed, which is the mean 
time for a Brownian particle to escape a confined domain through a small window p^[l6l.[T7j. However, 
all these computations rely on the assumption that the reaction volume is quite homogenous and has a 
shape close to a convex domain (no bottle neck). 

In photoreceptor outer-segment, a diffusing molecule needs to find a specific target site in a degenerated 
domain, where one dimensional length is much smaller than the others, and thus previous related to the 
small hole formula do not apply. This problem contains two difficulties: first, the target site occupies 
only a tiny portion of the boundary, and second, the diffusion occurs in a narrow domain. 

We shall now be specific and explain what is our goal in the context of phototransduction: Rod 
photor ecep tors are highly specialized biological devices that can detect a single photon absorption [isl 
[l9l. [20l. [2ll| . The photon absorption activates a cascade of chemical reactions in the outer segment, which 
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ultimately hyperpolarizes the cell [2l], [22, [23, [2J, [2^, [20| . The inner structure of a rod outer segment 
is very specific and can be considered as a cylinder that contains a densely packed stack of parallel 
and uniformly distributed discs (see Fig. [1]). The discs divide the outer segment into almost separate 
compartments that are loosely connected through a narr ow g ap between the disc perimeters and the outer 
segment membrane, which we refer to as the outer shell [27]. Compartments are also linked through disc 
incisures, however, since their impact is small 28] we will neglect them in first approximation . The 
chemical reactions involved in the early steps of phototransduction occur on the surface of the internal 
discs, and result in the activation of the phosphodiesterase molecule (PDE) via a G-protein coupled 
activation cascade [2l|, [13, [H, [2^. A photon-excited rhodopsin activates many transducin molecules, 
which bind to and thereby activate PDE. The number of activated PDE molecules following a single 
photon absorption was studied both experimentally and theoretically [l^, [s^, [HI, [13, [H, [s^l ■ We refer 
to PDE molecules that become activated via the phototransduction cascade as light-activated PDE. In 
addition to the transduction pathway, PDE can also spontaneously activate, leading to a non-vanishing 
background activity even in darkness [13, [1H| ■ 

Cytoplasmic diffusible cGMP molecules controlling the opening of cationic channels in the plasma 
membrane are hydrolyzed by activated PDE, and the reduction in the cGMP concentration leads to 
channel closure and photoreceptor hyperpolarization. From another chemical pathway catalyzed by 
guanyl cyclase (GC), a molecule attached to the disc surfaces and the outer segment membrane, cGMP 
molecules are synthesized from cytoplasmic GTP, a reaction which is calcium dependent. The magnitude 
of the photoresponse signal depends significantly on the number of closed ionic channels, and therefore 
on the drop in the cGMP concentration, which is controlled in part by the number of activated PDE and 
the rate of GMP hydrolysis of a single activated PDE. 

cGMP hydrolysis is characterized by two rate constants (id and (3subi which are both derived from 
experimental measurements [H, [13, [13, [13, [13, [H, [13 ■ Our goal here is to derive these constants from 
molecular considerations and biophysical theory, and thus obtain explicit analytical expressions. To 
understand at an intuitive level how these rates are defined, we recall that in most photoresponse models 
the cGMP concentration in the outer segment is well-stirred, a simplification that neglects diffusion and 
the complex geometry of the outer segment. The effective differential equation for the well-stirred cGMP 
concentration C{t) is [13, [13, [13| 

jC{t) = a{t) - (3dC{t) - (3subPi* {t)C{t) , (1) 

where is the number of light- activated PDE molecules and a{t) the rate of cGMP synthesis. The 

term f3dC{t) accounts for cGMP hydrolysis due to spontaneous PDE activation, and (3subPi{t)C{t) due to 
light-activated PDE. Eq.[l]shows an important difference in modeling cGMP hydrolysis by spontaneously- 
and light-activated PDE: whereas (3d is the rate constant for the change in the cGMP concentration due 
to all spontaneously activated PDE in the outer segment, (3 sub denotes the change in the well stirred 
cGMP concentration due to a single light- activated PDE. 

In the literature, Pd and f3sub are considered as independent parameters, a hypothesis that is strength- 
ened by the finding that the experimental values for (3d and (3sub are around ls~^ and 10~'^s~^ respectively, 
and therefore are extremely different in appearance [13, [13, HE [13| ■ 

Based on the diffusional encounter process between a cGMP and an activated PDE molecule in the 
complex rod outer segment geometry, we obtain explicit estimates for (3d and (3 sub- Our analysis is 
motivated by several known results: First, in darkness, in average around one spontaneously activated 
PDE molecule is present in a single compartment [13, [131 , which suggests that diffusion is rate limiting 
for hydrolysis. Second, experimental data [13,[13| indicate that activated PDE is a nearly perfect effector 
enzyme and hydrolyzes cGMP with a very high efficiency, which also hints that cGMP-hydrolysis is 
limited by diffusion. Third, a diffusion limited hydrolysis reaction couples the cytosolic cGMP level most 
strongly to the activation status of PDE, which is at the basis of photoreceptor adaptation [13, [13|- 

One of the main results of this paper is formula [5^ 

27rpA^+ 8 

which relates (3d to the spontaneous PDE activation and deactivation rates /x+ and the PDE surface 
density p, the effective reaction radius a, the radius R of the outer segment, and the cytoplasmic cGMP 
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diffusion constant Dcg- Furthermore, by comparing this purely diffusional cGMP hydrolysis rate to 
experimental results, we can estimate the impact of the details of the chemical hydrolysis reaction. 

The paper is organized as follows: we first determine the rate constant of cGMP hydrolysis due to a 
single activated PDE as a function of the cGMP concentration, the cGMP diffusion constant and the 
geometrical structure of the outer segment. Using this result, we then compute the analytical expressions 
for j3d and fisub- We find that (3d is proportional to the mean number of spontaneously active PDE 
molecules in a compartment and not in the whole outer segment. We compare our analytical estimations 
with experimental measurements, and find good agreement. Our analysis suggests that the main reason 
for the discrepancy between [3d and (3sub is their incompatible definitions. By deriving j3d and Psub from 
molecular events, we show that they are no longer two independent parameters. Finally, we use our 
analysis to model the spatio-temporal time course of a photoresponse in a transversally well-stirred outer 
segment. 



II. RATE OF CGMP HYDROLYSIS BY ACTIVATED PDE 



In this section, we estimate cGMP hydrolysis rate constant by a driven by single activated PDE 
molecule P* when diffusion is the limiting step. Later on, we use this result to derive expressions for j3d 
and (3 sub- To illustrate our approach, we start with the molecular model for cGMP hydrolysis: 

cGMP + P* ^ cGMP-P* P*+GMP. (3) 

A cGMP molecule binds to a P* molecule with a forward rate kf and forms an intermediate complex 
cGMP-P* . This complex can either dissociate with a backward rate fcf,, or cGMP becomes hydrolyzed to 
GMP with a rate k2- We are interested in the rate kh by which cGMP molecules are hydrolyzed, which, 
at steady state, balances the cGMP production rate. In the restricted rod outer segment, the overall 
forward binding rate is kfGc, where Gc is the number of cGMP molecules in a single compartment. As an 
example, in darkness, Gc is roughly in the range 100-1000, depending on the radius of the outer segment 
[2^ 1. From Eq. [3l using the overall forward binding rate and Michaelis-Menton approximation, we obtain 

k2kfGcP* ... 

kh = , ^ , " , ^ • (4) 
fc2 + fcfo + kfGc 

In the physiological range of cGMP concentrations, we assume that ^2 ^ kfGc, which implies that the 
hydrolysis of cGMP-P* proceeds much faster compared to the formation of a new complex. Furthermore, 
since P* hydrolyzes cGMP with very high efficiency @, H^), we suppose that k2 ^ fc?,, and therefore 
neglect the backward rate. Under these circumstances, Eq. [4] reduces to 

ku = kfGcP* , (5) 

which has exactly the form of the hydrolysis term in Eq. [TJ Eq. [S] can be formally obtained by setting 
k2 — 00, which means that cGMP hydrolysis occurs instantaneously after the formation of the the complex 
cGMP ■ P*. In contrast, if we assume that k2 is small (^2 ^ kfGc and kt, <C kfGc), using Eq. [4l this 
implies that hydrolysis proceeds independently of the cGMP concentration with a rate k2P*, a scenario 
that is not experimentally supported [23j . 

For large values ^2, the cGMP hydrolysis rate in Eq. Ois determined by the forward binding rate kf, 
whose value depends on two parameters: the encounter rate ke of cGMP molecules with the P* site, 
and the probability p that cGMP ■ P* is formed upon encounter. The probability p depends on (largely 
unknown) molecular properties of cGMP and activated PDE. In order to extract the impact of diffusion 
on cGMP hydrolysis, we set p = I and presume that the complex cGMP ■ P* is formed each time a 
cGMP molecule encounters P* . Finally, by neglecting the molecular details of activated PDE, we do 
not distinguish between spontaneously- and light-activated PDE, and we consider only activated versus 
non-activated PDE. If mainly diffusional issues are relevant for cGMP hydrolysis, then the catalytic 
activities of spontaneo usly - and light-activated PDE should be very similar, as was already suggested by 
experimental findings [30] ■ 
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Symbol 


Description 


L 


Length of a rod outer segment 


R 


Radius of a disc 


d 


Gap between disc and outer segment membrane 


I 


Distance between two adjacent disc 


Id 


Width of a disc 


ap 


Radius of a PDE molecule 




Radius of a cGMP molecule 


ap + acG 


Sum of the radii of a cGMP and PDE molecule 


P 


PDE surface density 


M+ 


Spontaneous PDE activation rate 


P- 


Spontaneous PDE deactivation rate 



TABLE L Description of the parameters used in the model. 



Before starting the analysis, we give range values for the main parameters: cGMP diffuses in the 
cytosolic volume Vcyto of the outer segment with a diffusion coefficient Dcg ~ 100/zm^/s~^ [ssL IstL IsoL l40l| . 
In contrast, PDE molecules are attached to the disc surfaces, where they diffuse with a diffusion coefficient 
DpDE ~ 0.8/iTO^/s~^ [2§|. The exact geometrical dimensions of a rod outer segment varies between 
species [23|, i41|] : for example, the length L and radius i? + d of the outer segment in a toad rod are 60fim 
and 3/xm, whereas in a mouse rod they are 20^m resp. 3/xm [1^. The longitudinal distance I between 
two adjacent discs (the height of a compartment) and the width of a disc Id (see Fig. [T]) vary around 
Ibnm [41]- The width d of the outer shell is comparable to / [4l|. The total number of compartments 
Nc = jT]- in the outer segment is of the order N^, ^ 10^. Finally, we assume that the radii ap of a PDE 
molecule and GcG of a cGMP molecule are both comparable to the radius of a rhodopsin molecule, which 
is around l-2nm [29| . The parameters are summarized in Table ID 



A 

! Z 




FIG. 1: Section through a cylindrical rod outer segment, containing a densely packed stack of parallel and 
uniformly distributed discs. The volume delimited by two adjacent discs is called a compartment. 
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A. Analysis of cGMP hydrolysis due to a single activated PDE 



To describe the time course of cGMP concentration in the outer segment, we consider different players: 
cGMP molecules are independent and diffuse freely inside the outer segment domain H.. Whenever a 
cGMP molecule hits the boundary area dflh occupied by the P* molecules, it becomes instantaneously 
hydrolyzed. The synthesis of cGMP occurs on the surface dfl — dflh with a rate aa{x, t). We account for 
these interactions by using the density C{x,t) of cGMP molecules at position x and time t, it satisfies 
the diffusion equation with the appropriate boundary condition 

DcGACix,t), for X en (6) 

—aa{x, t) , for a; e dn — dVlh (7) 
, for a; e dVLh . (8) 

We shall now study Eq. [S]for a single compartment ilc- 



D.alcix,t) = 
C(x,t) = 



Approximation of cGMP hydrolysis in a single compartment 




FIG. 2: (a) Elementary cylindrical compartment with a P* molecule located centrally on the upper surface, (b) 
Cross-section view of the compartment. A cGMP molecule is hydrolyzed when reaching the boundary of the 
volume Qa a^t r — a. 



We now consider a compartment f2c in which a single PDE molecule is activated on either one of the 
two disc surfaces (in FigS^ P* is attached to the upper surface). In our approximation, cGMP hydrolysis 
rate is given by the flux Jh of cGMP to the surface area dflp* occupied by a P* molecule. To compute 
Jh, we shall make some approximations: 

We consider a uniform and time independent cGMP synthesis rate a^- Because cGMP synthesis is 
calcium dependent, this corresponds to a calcium clamped outer segment or at equilibrium ( this is the 
case in darkness). Furthermore, because the height I of a compartment is around a few nm, and much 
smaller compared to the radius R ^ /im, the time scale for longitudinal equilibration I^/Dcg is much 
shorter than the one for radial equilibration E? /Dcg- Hence, newly synthesized cGMP molecules 
at the surface quickly equilibrate in longitudinal direction before encountering a P* molecule, which is 
usually located far away compared to the compartment height / (except for the negligible amount cGMP 
synthesized in direct neighborhood of P*). This scenario is equivalent to having cGMP synthesized inside 
the compartment, and we therefore replace cGMP synthesis on the surface by synthesis inside the volume. 
The volume synthesis rate a.^ is linked to by 



(9) 
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where the factor 2 accounts for the two disc surfaces enclosing a compartment. With a volume synthesis 
rate ay, the diffusion equation for cGMP is 



DcG^C{x, t) + , for X E ^c- 



(10) 



The boundary conditions are given by Eq. [8]and Eq. [7]with aa{x,t) — 0. 

Because cGMP diffuses much faster than PDE (Dpde ^ Dcg) we neglect PDE motion 0, ll^, IS- 
Since at leading order approximation the exact position of the activated PDE is not relevant [44, 45, 4q|, 
we position P* at the center of the disc. We will discuss this issue in more detail in section IIIII In 
addition, we approximate cGMP molecules by infinitesimal points, and use the effective reaction radius 
a = ap + flcG for a P* molecule pol. IHIisIl . 

Because the effective diameter 2o ~ 6nm of the boundary area dflp* occupied by a P* molecule is 
comparable to the compartment height I 15nm, and the radius R ~ Sfim is much larger than a and 
I, the main limiting factor for cGMP hydrolysis rate is the speed by which cGMP molecules find P* . 
We note that once a cGMP molecule enters into a neighborhood of dflp*, since 2a is comparable to I, 
it has a high probability to hit d^lp* and become hydrolyzed. In a first approximation, we model the 
hydrolysis reaction by assuming that a cGMP molecule entering the small cylindrical volume Qa (given 
in cylindrical coordinates by r < a) above or below dflp-^ is instantaneously hydrolyzed (see Fig. [5]). The 
corresponding boundary condition is 



C{x, i) = , for r = 



(11) 



This condition leads to an overestimation of the true hydrolysis rate because cGMP molecules entering 
the domain Qa can as well leave this region without touching the surface 9ilp. . However, in appendix 
lAl we show that the overestimation is in the range of a factor 2 (see also the discussion in section ITTl|) . 

Having discussed the approximations, we shall now proceed to estimate the cGMP flux Jh into dfta- 
We are particularly interested in Jh as a function of the cGMP concentration inside the compartment. 
From there, we will extract cGMP hydrolysis rate due to a single activated PDE. With this result, we 
will then derive the expression for /J^; (see Eq. [Ij . Using the cylindrical symmetry, Eq. [TO] reduces to 



|C(M) = Dj-§-r§-Cir,t)+ay. 
C{r, t) = for r = a . 



(12) 
(13) 



Integrating Eq. [12] over the compartment volume yields an equation for the time dependent number of 
cGMP molecules Gc{t) in Qc, 



dt 



Gc{t) = -JrH) - Jh{t) + a^TTR^l 1 



a 



(14) 



where 



Gcit) = 2Trl / C{r,t)rdr, 

J a 

Jh{t) = 2TialD,G^C{r,t) 



Jpit) = ^2nRlD,G^Cir,t) 
or 



r=R 



(15) 
(16) 
(17) 



The fiux JR{t) is maintained by cGMP molecules that diffuse between compartments. To derive an 
expression for J^, we consider the steady state regime where the flux Jp is given. The steady state 
concentration C (r) obtained from Eq. [T^] is given by 



C{r) 



ayirR^ 
2ttDcg 



In 



2i?2 



Jr 
2tiDcgI 



■In 



(18) 
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To obtain the number of Gc molecules inside a compartment, we insert Eq.[TH]into Eg. llSl and for ^ 1, 
we obtain : 

o, i?^ /I. fR\ 3\ . fl, fR\ 1 



We define the times ti and T2 and the corresponding rates fci and ^2 as 

n = 7^ = l^(^ln(-)-|V (20) 



1 
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we can rewrite expression [19] as 

G. = ^^-:^. (22) 

fci fc2 

At steady state, the value of Jh is fixed by the balance of fluxes, and Eq. [T3] gives for ^ 1 

Jh = -JR+avirR^l. (23) 
Using Eq. [22]we can express as a function Gc and Ju, 

fci 

a„7ri?2? = kiGc + -t-Jr . (24) 

fc2 

Finally, inserting Eq. [24] into Eq. [23] yields 

A = fciG, + -^i^Jfl. (25) 

Formula [23 gives the steady state hydrolysis rate Jh as a function of Gc and Jr. This result depends 
strongly on the diffusional and geometrical properties of the microdomain. Whereas Eq. [53] gives a direct 
expression for Jh as a function of the synthesis rate av, and does not involve diffusion, Eq. [^^[is related 
to ay indirectly via the value of Gc, and therefore involves diffusion. 

In appendix [X] we obtain an interpretation for the two times ti and T2, and, thus, for the rates fci and 
fc2: Ti (see Eg. IA2l| ) is the mean time for uniformly distributed cGMP molecules to reach the absorbing 
boundary at r = a, given reflecting boundary conditions at r = i?; T2 (see Eq. lA17] for r = i?) is the mean 
time to reach r = a, when the initial position is uniformly distributed at r = i?. In reality, there is no 
reflecting boundary at r = _R, however, a vanishing flux Jr is mathematically equivalent to a reflecting 
boundary condition at r = _R. 



B. Derivation of the rate constant Pd for spontaneous PDE activation 

In darkness, spontaneous PDE activation leads to a uniform cGMP hydrolysis in the outer segment 
[23 , [25 , |30| with an overall hydrolysis rate (see Eq. [T] integrated over the cytoplasmic volume) 

Jd,os — PdGos , (26) 

where Gos is the total number of cGMP molecules in the outer segment. To derive an analytical expression 
for the rate constant /3d, we start from Eq. [^S] Because spontaneous PDE activation occurs uniformly 
throughout the outer segment, apart from fluctuations, the flux Jr between compartments vanishes in 
darkness. Thus, the steady state hydrolysis rate Jh of a single P* molecule given in Eq.[5S]can be written 

as 



Jh = fciGc 



(27) 
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To obtain the dark hydrolysis rate Jd^c per compartment, we have to further consider the mean number of 
spontaneously activated PDE molecules P* ^ in a compartment. As long as the number P* ^ is small and 
the P* molecules are geometrically well separated [47], the rate Jd,c increases linearly with P* ^. Hence, 
we obtain 



Jd,c = fciP*cGc. (28) 

The hydrolysis rate in the whole outer segment Jd,os is obtained by summing Jd^c over all Nc compart- 
ments. Since Gos — NcGc well approximates the total number of cGMP molecules in the outer segment 
(the volume inRdL of the outer shell is negligible compared to the volume ttR^INc of all compartments), 
we obtain 



Jd,os = kiPlcGos ■ (29) 
Finally, by comparing Eq. [21] with Eq. [55] and by using Eq. [501 we obtain 

^'^-^^^^=-^Ii^|F3^-- 

We conclude that Pd is determined b y th e mean number of spontaneously activated PDE in a compart- 
ment, and not in the outer segment [30|. Furthermore, we will now relate Pd to the spontaneous PDE 
activation rate , the deactivation rate fj,- , and the PDE surface density p. The number of PDE on the 
disc surfaces attached to a single compartment is Pc = 2ttR^p, and P* ^ is given by 

p;^ = = 27ri?2p— • (31) 

fi- p- 

Together with Eq. [20] and Eq. [30] we obtain the final expression 

P,^D,^^^ 1 . (32) 

p^ 41n(f)-3 ^ ' 



C. Effective set of equations to model cGMP dynamics 



By generalizing our previous results, we shall now derive an effective set of equations to model cGMP 
dynamics following a photon absorption. Since a photon absorbtion transiently generates an elevated 
amount of P* molecules inside the affected compartment, it induces an increased cGMP hydrolysis and 
a cGMP gradient in the outer segment. In this case, the fluxes Jr between compartments are no longer 
zero after a photon absorption. 

We start the derivation by extending the equilibrium expression for J/,, given in Eq.[25]to time dependent 
situations. Because free cGMP diffusion is fast, cGMP equilibrates quickly inside a compartment. In 
contrast, Gc{t) and JR{t) fluctuations are determined by the effective longitudinal diffusion between the 
compartments, which is strongly hindered by the compartmentaHzation of the outer segment [27l. I40ll48j . 
Thus, we consider that Gc{t) and JR{t) fluctuate on a slower time scale compared to the equilibration 
time scale inside a compartment. Under this condition, a flrst approximation of the time dependent 
hydrolysis rate Jh (t) of a single P* molecule is given by the equilibrium expression in Eq. [25] with time 
dependent Ji?(i) and Gc{t): 

Ju{t) = fciGe(t) + ^l^Mt) . (33) 

We note that the expression for Jh given in Eq. [55] can be extended to time dependent cases, whereas 
this is not possible starting from Eq. [53] 

To obtain the set of equations for the time dependent number of cGMP molecules inside a compartment, 
we start when a photon is absorbed in compartment tiq, while the other compartments n — 1 . . . Nc, 
n ^ no remain unperturbed. In the regime considered here, cGMP hydrolysis depends linearly on Pi*{t). 
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Using Eq. [TU the equation for the number ci""* (t) of cGMP molecules in a compartment n is given by 
{Sn.no is Kronecker-Delta) 

|G(")(t) = - Jl,"'(t) - (P;,, + Pnt)Sn,no)jt\t) + a^nRH . (34) 

Inserting the expression for jj^"^ (t) given in Eq. 1331 and using the definition of f3d in Eq. 1301 we obtain 

^Gi-Ht) = - (l + (PL + Ptmn,n.)^-^) 4'\t) 

-fidGi^'Ht) - hPi*{t)Sn,n„Gi"'\t) + a^TiRH . (35) 

By approximating the transversally well stirred cGMP concentration in a compartment by C''^\t) « 
'^ttbH ' using Pick's law, the fluxes J^ji\t) are approximated by 

(„) ^^^ G("+i)(t)-G(")(t) , C("-i)(^)-C(")(0 
it) = -7^,G27ri?d(^ + 

^ 2d ( G^J'+'\t)~Gi"\t) , Gi''-'\t)~G^^\t) \ 

= """^^^a 1, TTTd + TTL J ■ 

Egs. [551 and [5^ constitute a close system of equations for the Gi"\t) (which can be transformed into 
equations for the concentrations C^'^\t)). Furthermore, Eq. [351 models the impact of spontaneously- and 
light-activated PDE in an equivalent way. The simulation in Fig. [3l shows the time course of the number 
of cGMP molecules, scaled with respect to the dark equilibrium value, after the absorption of a photon 
at time t = in the middle of the outer segment. The parameters for the simulation are suitable for a 
toad rod [11] . The input function Pj* (t) is obtained using the set of equations published in [s^] . 



D. Derivation of the rate constant /3sub in a well stirred outer segment 

We now derive an analytic expression for the rate constant Psub using the approximation of a well 
stirred outer segment [l^. Since the volume of the outer shell is negligible compared to the combined 
volume of all compartments, the total number of cGMP molecules in the outer segment is 

Gosit) = Y.Gi-Ht). (37) 

n=l 

By summing Eq. [551 over all compartments, and using that J2n=i '^r\^) ~ ^ neglect the cGMP 
molecules in the outer shell), we obtain 

j^Gosit) = -^^j^Pi*{t)/^"\t) - PdGosit) - fciP;(OG("«)(i) + a.TTR^lNe . (38) 

In a well stirred outer segment we have Gi"'\t) = Gos{t)/Nc. By further neglecting the terr 

1 — 

k2 



''^^^Pt{t)J^f^"\t) (the flux Jr vanishes in a well stirred outer segment), we get 



^Gos{t) = -PdGosit) - ^Pt{t)Gos{t)'^a^^RHN,. (39) 
at Nc 

Finally, dividing Eq. [39] with the cytosolic volume Vcyto ~ ttR^INc yields the standard equation for the 
well stirred cGMP concentration C{t) — Gos{t) /Vcyto, 

jC{t) = -ddGlt) - ^Pt(t)C{t) + a„ . (40) 
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FIG. 3: cGMP dynamics after a photon absorption at time t = in compartment n — 1000 {Nc = 2000). The 
simulation is performed using Egs. 1351 and 1361 The cGMP concentration is scaled with the equilibrium value, (a) 
Time dependent cGMP level averaged over the outer segment, (b) cGMP level per compartment for various time 
points, (c) Number of light-activated PDE molecules obtained using the equations published in [3^ . 



By comparing Eq. [TD] with Eq. [T] we obtain for Psub the expression 

fci Pd 



(41) 



Since Nc is of the order 10'^, it follows that f3sub is much smaller than f3d- Using Eq. [^D] for ki and 
Vcyto ~ nR'^lNc, Eq. [lT]can be written as 



sub 



Vcyto 41n(f) -3 



(42) 



By comparing exDressioi il42l with the standard definition of Psub given by [22, (we neglect cytoplasmic 
buffering for cGMP ^,^) 



^sub 



Km^AvVcyto 



(43) 



11 



we obtain a new formula for given by (Nav is the Avogadro number) 



ksub _ NAyVcytpki _ 8tt DcgINav , . 



III. COMPARISON WITH EXPERIMENTAL RESULTS 



To validate our computations, we now compare our analytical results for and /?s„f, (Eq. [30] and 
Eq. [4T|) with experimental measurements [l^, [23, [sO, [l^l- We start with Pd- Using data available for 
toad rods, p = lOOum^^, R = 3fim, /i+ = 4 x 10~^s~-^, ^_ = 1.8s~^, Dcg = lOO^m^/s^^, a = 3nm, 
Z = 15nm [i^ [13, 113], and inserting these values into Eq. [20l Eq . [3T] and Eq. [30l we obtain ki = 3.6s~^, 
P;^ = 1.26 and 



I3d = A:iP*, « 4.5S-1 . (45) 

This analytic result has to be compared to the experimentally found value (3d ~ ls~^ [2^ , which is 
approximately four times smaller than this prediction. Eq. [35] shows that (3d depends only logarithmically 
on the compartment radius R, and thus it is very similar across species that differ mostly on the radius of 
the outer segment, in agreement with experimental findings [23j . The discrepancy between our theoretical 
prediction and the experimental value for (3d can be attributed to several factors: 

1. We made the assumption that a cGMP molecule already becomes hydrolyzed when reaching the 
inner cylinder fia at r — a. Thus, the time ti in Eq. [20] is shorter than the true time needed to 
arrive at the P* site. Hence, Eq. [20] overestimates the hydrolysis rate. In appendix [X] we derive 
an accurate estimate for the mean time t a cGMP molecule reaches the P* site located on the 
surface of a compartment (see Eq. IA23p . Compared to n (Eq. [20ll . the new estimate for r includes 
specifically the mean time Tq a cGMP molecule starting on the boundary of ^la reaches the P* 
molecule on the surface. By considering the additional time Ta, we replace ti and T2 with the more 
accurate expressions fi — Ta + ti and f2 = Tq + T2 . Accordingly, the rates ki and k2 have to be 
replaced by ki and k2, given by 

k = ^ = ^—, ~k2 = ^ = ^—. (46) 

n Ta + Ti T2 Ta + T2 

For toad rod values with l/a ^ 5 and R/a ^ 1000, and by using Eq. IA23l with g(5) « 2.9 (the value 
g{5) is obtained from Fig. [51d), we find that ki « 0.5/ci. By using ki instead of fci in Eq . [30l we 
obtain the new estimation 

(3d = hP:,c = 2.25S-1 , (47) 
which is closer to the experimental observation. 

2. Our assumption that every encounter between cGMP and P* results in cGMP hydrolysis will cer- 
tainly lead to an overestimation of the hydrolysis rate. Moreover, since we neglected the molecular 
details of the hydrolysis reaction, this will also induce an error. Nevertheless, since our analytic 
result for (3d is very close to the experimental finding, we conclude that cGMP hydrolysis by P* 
has to be largely diffusion limited, and in addition has to be quite efficient, such that nearly every 
encounter between cGMP and activated PDF leads to a hydrolysis reaction. This is supported 
by the experimental observations that activated PDF hydrolyzes cGMP with very high efficiency 

3. Uncertainties in the experimental values for DcGi /^+ ^-nd involved in the computation of 
(3d, introduce ambiguities in our analytical prediction. For exainple, there is still considerable 
disagreements about the exact value of the diffusion constant DcG [S HSj IHj ISI ■ Furthermore, at 
first approximation, we used for the effective reaction radius a the sum of the molecular radii of a 
PDF and cGMP molecule. A more precise value for a will affect Pd in Eq. [47] mainly via Ta, since 
Ti depends only logarithmically on a (see Eq. [32] 
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4. The value of jSd was computed by fixing the position of P* at the disk center and neglecting 
possible fluxes between compartments. In general, spontaneous PDE activation and diffusion leads 
to P* positions that are uniformly distributed over the disk surface, and different P* positions in 
neighboring compartments induce small fluxes. We left open here the computation of the variance 
of the cGMP hydrolysis rate constant coming from random locations of P* molecules. However, the 
P* position should not much influence the rate constant for cGMP hydrolysis: The rate constant 
is determined by the MFPT of a cGMP molecule to find the P* target. Outside a small boundary 
layer around P* (the radius of the boundary layer is of the order of the reaction radius a), the 
leading order term of the MFPT in dimension 2 depends only logarithmically on the distance 
between cGMP and P* , and in dimension 3 it is a constant [H, . Hence, since almost all cGMP 
molecules are outside the boundary layer, the exact position of P* is not important for their mean 
time to hydrolysis. We conclude that our expression for (3d should remain a valid approximation at 
first order, even when considering random P* positions. 

We shall now compare expressions Eqs. I41I42I for Psub, and the ratio (Eq. with experimental 
measurements. Eg . |411 reveals that /3sub is a factor NcP* smaller than Pd- Since Nc is of order 10^ and 
P* ^ of order 1 — 10, this agrees with the experimental findings that f3sub is around 10^ — 10^ times smaller 
than Pd [H, [2^ . From Eq. [331 we obtain the prediction 

Km 

This estimation can be further improved by using the rate fci = 0.5fci instead of ki, giving 

^^"^ - 4.6 X loHr^s-^ 



which has to be compared to « 2.2 x 10^M~^s~^ obtained from experiment [s^l- It is important 
to note that our analytic results for and Pd (see Eq. [37]) are both around two times larger than the 
experimental findings, which is an indirect confirmation of our assumption that f3d and f3sub (note that 
ksub/ Km is proportional to Psub) are not two independent rate constants, but can be derived from the 
same underlying hydrolysis reaction. Despite of the encouraging results, we would also like to indicate 
some difficulties related to the definition and derivation of the parameters Psub and First, we 

extracted the formula for using the expression for Psub given in Eq. [321 This approach is problematic 
because the definition of Psub involves the assumption of a well stirred cGMP concentration during a 
photoresponse, which is not very accurate (see Fig. [3] and [13, [s^])- Second, if diffusion limits the rate of 
cGMP hydrolysis in the physiological range, the experimentally observed value for does not reflect 
an intrinsic property of the chemical reaction. Instead, it depends strongly on diffusional and geometrical 
details, and, therefore, on the experimental setup. For example, measurements of the Michaelis constant 
Km were performed using fragments of disrupted rod outer segments with a length only a fraction of 
the intact outer segment length [H, . Eq. [40] shows that the rate of cGMP hydrolysis increases with 
decreasing fragment length L (since Nc ~ L). Thus, the apparent value of the Michaelis constant Km 
{ksub is assumed to be a true constant) that is needed to fit the rate of cGMP hydrolysis will be higher in 
a suspension containing large fragments compared to a suspension with small fragments, as it has been 
observed [3^.|49|. 

In this work we have assumed that cGMP hydrolysis in the physiological range is diffusion limited, 
and is independent of whether PDE is spontaneously- or light-activated. The agreement between our 
theoretical results and experimental measurements indicates that the large disparity between psub and f3d 
is largely due to their definition, and not due to biochemical differences. For example, in the effect of 
spontaneously- and light- activated PDE was modeled using two very different rates k — 0.042uM" 



and k* — 110/iAf ^^s^^. We will show now that the large discrepancy between k and k* in [27[ essen- 
tially originates from modeling needs. Indeed, cGMP hydrolysis by spontaneously activated PDE was 
modeled as k[PDE]c[cGMP], where [PDE]^ is the surface concentration of PDE. In contrast, hydrolysis 
by light- activated PDE was modeled as k*[PDE*]cr[cGMP], with [PDE*]c as the surface concentration 
of light-activated PDE. By introducing the mean surface concentration of spontaneously activated PDE, 
[PDE*]„ = [PDE]„i^, we rewrite k[PDE]„[cGMP] as k^[PDE*]„[cGMP], which now has the same 
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form as k*[PDE%[cGMP]. Inserting the values m+ = 4 x IQ-^s-^ and ^_ = l.Ss^^ found in (s^l, we 
obtain kj^ — 189/iAf^^s^^, which is now comparable to k* — 110/iM~^s~^. We conclude that model- 
ing cGMP hydrolysis by spontaneously- and light-activated PDE in a similar way involves comparable 
parameters, indicating that hydrolysis may be indeed independent of whether PDE is spontaneously- or 
light-activated. 

IV. CGMP HYDROLYSIS IN CONES 

After having discussed in detail cGMP hydrolysis in rods, we now briefly explore hydrolysis in cones. 
Similar to 48], our analysis for rods can be adapted to cones. Unlike rods, cones do not contain disc in 
the outer segment. However, the membrane invaginations in cones can be modeled similarly to discs in 
rods. Since the radius of the cone outer segment decreases from the bottom versus the top, we can adapt 
our formulas to cones by replacing the disc radius R with a compartment dependent radius i?„. Thus, in 
cones, the rates fci and ^2 depend on the compartment n. Therefore, the response to a photon absorption 
in cones varies on the location where the photon has been absorbed. Since kiP* ^ depends logarithmically 
on the compartment radius i?„ (see Eq. I^O)) , we suggest that the value for the dark hydrolysis rate fid in 
cones should be of the same magnitude as found in rods, see also [s^ . 

V. SUMMARY AND DISCUSSION 

In this paper, we have studied the rate constant of cGMP hydrolysis by activated PDE in rod and cone 
photoreceptors. Our analysis is based on the assumption that cGMP hydrolysis is diffusion limited and 
determined by the encounter rate between cGMP and activated PDE. We derived an explicit formula 
for the rate constant of cGMP hydrolysis by a single activated PDE molecule as a function of the 
confined outer segment geometry and the cGMP diffusion constant (Eq. [27]). Our calculation takes into 
account the complex structure of the rod outer segment, uniformly divided by a stack of parallel discs 
into homogenous microdomains, called compartments, and coupled to each other via cGMP diffusion. 
We obtained analytical expressions for the rate constants (3d and /3s«6- In addition, we give a set of 
effective equations that allow to model the transversally well stirred cGMP concentration after a photon 
absorption. 

Interestingly, we found that only the amount of spontaneously activated PDE in a single compartment 
is needed to calculate the dark hydrolysis rate Pd (see Eq. [50)1 . This result differs from where the 
compartmentalization was not considered, and all spontaneously activated PDE in the outer segment 
additively contribute to Pd- Because the number of spontaneously activated PDE in the outer segment 
is by a factor Nc ~ 10'^ larger compared to a single compartment, the catalytic activity of an excited 
PDE in [s^l was estimated much lower compared to what we found here. We computed the PDE activity 
(given by the rate fci in Eq. to be around ls~^, whereas in |33| it is around lO'^s""'^. Using the 
rates for spontaneous PDE activation and deactivation (30| . we estimate that the average number of 
spontaneously activated PDE molecules in a single compartment is around one. Together with our result 
for the PDE activity, this naturally explains the experimental value j3d ~ Is^^. For the derivation of (3d, 
it was essential to assume that cGMP hydrolysis occurs locally at the activated PDE site. In contrast, if 
cGMP hydrolysis occurred uniformly over the disc surface, then the experimental value for (3d could not 
be recovered without introducing additional adjusting parameters (see appendix [B|) . 

We have derived a set of equations (Eqs. [3S] and [55)) that allow to calculate the time course of the 
transversally well stirred cGMP concentration following a photon absorption. These equations model 
cGMP hydrolysis by spontaneously and light-activated PDE in a similar way. Under the assumption 
that cGMP concentration in the outer segment is well stirred, we derived an expression for the rate (3sub 
fEq. [iTB^ . and the ratio ^ (Eq.[ll|). 

Eq. [51] connects (3 sub to (3d and gives a direct explanation why (3 sub is found to be so much smaller than 
(3d- Our result suggests that the large discrepancy between (3d and (3 sub is largely due to their definitions: 
(3d incorporates the effect of all spontaneously activated PDE in the outer segment, while (3 sub accounts 
for only a single light- activated PDE. 
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APPENDIX 

APPENDIX A: MEAN TIME TO HYDROLYSIS IN A COMPARTMENT 




(b) 



FIG. 4: (a) Cylindrical compartment Qc with an activated PDE molecule located centrally on the upper surface, 
(b) The first time a cGMP molecule hits dQp* , when starting at position x in Qo, is given by the sum of two 
times. First, the time To for a molecule starting at x to arrive at the boundary dila, and second, the time Ta for 
the particle starting at dQa to arrive at dQp* . 



In this part of the appendix, we shall obtain a precise estimate for the mean time r a cGMP molecule 
starting uniformly distributed inside the cylindrical compartment flf. reaches the activated PDE molecule, 
defined as the small surface patch dflp* with radius r — a (see Fig 2]). The motion of the cGMP molecule 
is Brownian in flc. It is reflected all over the boundary except at dflp* , where is is absorbed. We denote 
by T^^) the random initial time a cGMP molecule starting at x d flc hits dflp*. Due to rotational 
mvariance, the Mean First Passage Time (MFPT) t{x) = E[r('') |a;(0) = a;] depends only on r and z. 
Using a cylindrical coordinate system x = (r, ip, z)), we decompose the domain fie into the inner cylinder 

D,a = {x e ^c\r < a} , (Al) 

and the hollow cylinder 

^o^n^-na^ {x\a <r <R}. (A2) 
We define the mean time r as the average over a uniform initial distribution in f2c, 

Using that fl^ ~ ilo + ^a, we can rewrite Eg. IA3I as 

_ f ,(,,w + ^M_^^ / r(rz)dV 

= [ ^(^^^W-MIt^ [ ^^r,z)dV--^, I T{r,z)dv], (A4) 



where i^pj- /^^ T(r, z)dV is the mean time for particles starting uniformly distributed in ilo to hit dilp* , 
and T(r,z)dV is the mean time to hit d^lp* for particles starting uniformly distributed in fla- 

Because particles originating from fio have to reach fia before hitting dH. p' , it is plausible (and can also 
be shown) that the mean time to dQ.p* for particles starting uniformly in is larger than the mean 
time for particles starting uniformly in f7a- Furthermore, for a <ti R, we have \Qa\ \^c\- From this, we 
finally obtain 

^ ^ T{r,z)dV -1- [ T{r,z)dV. (A5) 
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Thus, for a <C i?, the mean time t is well approximated by the mean time for particles starting uniformly 
in Qo to hit dflp* . 

We shall now estimate T(r, z), by considering the equation P. [soj 

/ (92 Id 92 \ 

^=G(^ + -^ + ^)r(r,z), = -1 0<z<;,0<r<i? (A6) 



, z — I ,r > a 





0, 


= -] 


T(r, 


^) 


= 0, 


oz 


^) 


= 0, 


oz 


^) 


= 0, 


^ ( 


^) 


= 0, 



= 0, r = R. 



We will first estimate the average time 



1 



T{r) = - j T{r, z)dz . (A7) 



Using the dimensionless variables 



x=-, y=-, T{x,y) ^ ^^T{r,z) , a = 2^0 = -, (A8) 



equation IA6I becomes 

dx dx dy 



I d d 92 

-^2;— + 1 ^(x,?/) = -a^ 0<?;</3,0<x<a;a (A9) 



92 \ 


= -a' 


f(x,?/) 


= 0, 


-TT'^yx.y) 
ay 


= 0, 




= 0, 


= 0, 



and Eq. |A7l 



1 



/3 



t{x) = - / T{x,y)dy. (AlO) 



We integrate Eq. IA9I over the variable y to derive an equation for t{x) for a; > 1. Taking into account 
the boundary conditions at y and y = (3, we obtain 

19 d 

-—x—t{x) = -o? , x>\ (All) 
X ax ox 



The solution is given 



with 



d ^ 

—t{x) = for x ~ Xa 
Ox 



r{x) = /(a, /3) + ^ ln(x) - \a\x^ - 1) , (A12) 



/(a,/3) = f(l) = i^%(l,2/)dy. (A13) 
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Hence, for r > a, we have 



r(r) 



i?2 



where we defined 



DcG 
Ta + To{r) , 

i?2 



To{r) 



i?2 



cG 



-/(«,/?), 



D. 



cG 



-In - 



4 i?2 



(A14) 
(A15) 

(A16) 
(A17) 



Eg. IA14] has an intuitive interpretation: the mean time r(r) for a cGMP molecule, uniformly distributed 
at r > a, is the sum of the mean time To(r) to the boundary r — a plus the mean time Ta from the surface 
dila to go to 9rip» (see FigH]). 

By averaging over a uniform initial distribution p — ^(^]^i_ai) ^o, the overall mean time r in Eg. IA5I 
is given by 



1 

D. 



"(r, z)dV = 2'Kp I T{r)rdr 



cG 



'4 In (a) - 3 + 4^2 



8-D, 



cG 



where we defined To as 



To = 2'Kp / To{r)rdr 



i?2 



1 - a2 



-4 In (a) + -a^ 



cG 



1 -a^ 



The leading order expansion of Tq for a ^ 1 is 

i?2 



I? 



cG 



1 



i? 



For a ^ 1 (see also Fig.lSJj), we can approximate /(a,/3) by /(0,/3) = g{l3) and obtain (/3 = ^) 

i?2 



i:'cG 



5 



a < 1 



Altogether, for a <IC 1, the mean time r in Eg. IA5l is given by 



D. 



cG 



9 I - 

a 



1 



R 



;in - 



(A18) 
(A19) 

(A20) 
(A21) 

(A22) 
(A23) 



To derive an explicit expression for /(a, P) and g{P) is a difficult mathematical problem. Nonetheless, 
we shall obtain some asymptotic limits for f{a,f3). For /3 — *■ (corresponding to Z 0), we have Ta — + 0, 
and therefore f{a, 0) = 0. For /3 oo, the time Tq diverges to infinity, and /(a, f3) oo. For a = R, 
corresponding to a = 1, we have = JTT^^ from which it follows that 



0^ 



(A24) 



Finally, the small hole theory f44l| predicts that when I ^ R <^ a (which implies that a ~ /3), the mean 
time r is asymptotically given by 



V 



(3. 



(A25) 
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By comparing Eg. IA25I with Eq. IA23I (ln(Q;) can be neglected compared to /3 for a ~ we obtain the 
asymptotic 



g(/3)~-/3, /3»1. 



(A26) 



So far, we have only an asymptotic expansion for /3 ^ 1. To explore a much larger parameter space, we 
decided to run Brownian simulations (10000 cGMP molecules are initially uniformly distributed over the 
lateral surface of dila) to estimate and J {a, (3). The numerical results for f{a,f3) are summarized in 
Fig. El Fig. [5)3 shows that f{a,f3) can be well approximated by /(0,/3) = g{P) for a < 0.05. 




6 
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°0 0.01 0.02 0.03 0.04 0.05 

a 



(a) 



(b) 



6 
4 



2 



4 6 

(c) 



10 



FIG. 5: Numerical evaluation of the function f{a,P) {a = a/R and /3 = I /a). Each data point is obtained using 
the Brownian simulation of 10000 cGMP molecules, (a) f{a, (3) for different values a and /3. For a = 1 we have 

f{l,(3) = ^ (Eg. IA24)) . (b) Same data as in (a), restricted to small a. (c) Plot of /(j^,/?) ~ g{P) (same data 
as in (a) and (b)). The dashed curve represents the asymptotic j/S (Eg. IA26)l . achieved for a — > and /3 — > 00 . 
For values /3 < 10 used in the simulations, the behavior of g(/3) is close, but not yet in full agreement with ^/?. 
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APPENDIX B: MODEL WITH UNIFORM HYDROLYSIS ON THE DISC SURFACES 

In this section, we consider a model that is based on cGMP hydrolysis occurring uniformly on the disc 
surface (see also [27, 37]). This is very different from the situation presented within the main body of the 
paper, where we assumed that cGMP hydrolysis occurs locally at the P* site. We show now that uniform 
cGMP hydrolysis leads to a dark rate constant proportional to Dcg/P' ^ 10^s~^, which is very different 
from Eq. [30] (see also section Hill) . Thus, in order to account for the experimental value (3d ~ Is^^, one 
is forced to introduce a small adapting parameter k/j. 

We will analyze two different scenarios: In one situation, synthesis and hydrolysis of cGMP are both 
modelled by boundary source terms. In another, only hydrolysis is modelled by a boundary source term, 
while synthesis occurs uniformly within the cytoplasmic volume. 



1. Model with boundary source terms for hydrolysis and synthesis 



The reaction-diffusion equation for cGMP concentration C{z,r,t) inside a compartment is given by 

dt 

dC{z, r, t) 



-Ciz,r,t) = D,G^Ciz,r,t)^ 



DcG 



dz 

dC{z, r, t) 
dr 



5=0 
'.=1 



-KhC{z, r, t) ,=0 + a„{t) 

z — l 



= 



(Bl) 
(B2) 
(B3) 



The steady state expressions for the concentration C of Eq. lBlj and the hydrolysis rate Jh are (Gc = ttR^C) 

C ^ ^, (B4) 
Jh = 2TTR^KhC (B5) 



Because hydrolysis and synthesis are both modelled by fluxes originating from the same boundary, the 
equilibrium only reflects the balance of the fluxes, and does not involve cGMP diffusion. 



2. Model with boundary source term for hydrolysis, and a volume synthesis rate 

To include cGMP diffusion (see section III A[) , we now model cGMP synthesis by a uniform volume 
production rate av{t) = 2arj{t)/l. This model avoids the problems arising when cGMP synthesis and 
hydrolysis are both modelled by surface fluxes originating from the same boundary. The equation for the 
cGMP concentration reads 



d 

— C{z,r,t) = DcGAC{z,r,t) + ay{t), 
at 



-D,G 
DcG 



dC{z, r, t) 



dz 

dC{z,r, t) 



dr 



z = 



r=R 



-KhC{z,r, t) 



i = 

z=l 



(B6) 
(B7) 
(B8) 



The steady state solution of Eq. IB6I is 



(B9) 



where we introduced the parameter [3 and the concentration C as 



DcG 



'DcG 



(BIO) 
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At equilibrium, the number of cGMP molecules Gc and the hydrolysis rate Jh in a compartment are 







Gc = TiR^ / C{z)dz = C ( — + — ) ttR^I , (Bll) 



CcG / 1 , 1 ' ' 



P yi2 2/? 



Gc. (B12) 



Contrary to Eg. IB5[ the flux Jh in Eg. IB12I depends on cGMP diffusion constant. In the limit /3 ^ oo 
('t/i — > oo) (perfectly absorbing boundaries), we obtain Jh = Gc/t, where r — ^^-^ ^ is the mean time for 
a molecule to reach the boundaries at z = Z or z 0. On the other hand, in the limit — > 0, we obtain 
Jh = ^^ff^'2' PGc = 2^Gc- Thus, we recover the expression given in Eg. IB5I 

Eg. lB12l is formally equivalent to Eq. [27l however, the physical content is very different. Jh in Eg. IB12I 
is proportional to the rate by which cGMP molecules collide with the disc surfaces, which is of the order 
^ lO^s^^. In contrast, Jh in Eq.[27|is determined by the rate by which cGMP molecules find P* , 
given by Dcg/R^ ^ Is^^. In order to obtain Pd ~ Is"^ from Eq. IB121 one needs a small value for the 
adapting parameter Hh- 
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